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We report a new collective dynamical state, namely the mixed mode oscillation suppression state 
where different set of variables of a system of coupled oscillators show different types of oscillation 
suppression states. We identify two variants of it: The first one is a mixed mode death (MMD) 
state where a set of variables of a system of coupled oscillators show an oscillation death (OD) 
state, while the rest are in an amplitude death (AD) state under the identical parametric condition. 
In the second mixed death state (we refer it as the MNAD state) a nontrivial bistable AD and a 
monostable AD state appear simultaneously to different set of variables of a same system. We find 
these states in paradigmatic chaotic Lorenz system and Lorenz-like system under generic coupling 
schemes. We identify that while the reflection symmetry breaking is responsible for the MNAD 
state, the breaking of both the reflection and translational symmetries result in the MMD state. 
Using a rigorous bifurcation analysis we establish the occurrence of the MMD and MNAD states, 
and map their transition routes in parameter space. Moreover, we report the first experimental 
observation of the MMD and MNAD states that supports our theoretical results. We believe that 
this study will broaden our understanding of oscillation suppression states, subsequently, it may 
have applications in many real physical systems, such as laser and geomagnetic systems, whose 
mathematical models mimic the Lorenz system. 

PACS numbers: 05.45.Xt, 05.45.Gg 


I. INTRODUCTION 

The suppression of oscillation is an important emer¬ 
gent behavior shown by coupled oscillators, and it has 
been extensively studied in a variety of fields such as 
physics, biology, chemistry and engineering [l|. Ampli¬ 
tude death (AD) and oscillation death (OD) are the two 
distinct oscillation quenching processes. In the AD state, 
oscillations cease and the coupled oscillators reside in a 
common stable steady state; thus, it induces a stable ho¬ 
mogeneous steady state (HSS) which was unstable 
otherwise. On the other hand, OD is a much more com¬ 
plex and completely distinct phenomena than AD [l|. 
In OD, oscillators occupy different coupling dependent 
steady states and thus give rise to a stable inhomoge¬ 
neous steady state (IHSS). AD is a widely studied topic 
as a control mechanism to suppress undesirable oscilla¬ 
tions in laser application @|, neuronal systems |||, elec¬ 
tronic circuits [5j, etc. On the other hand, OD has strong 
connections and influence in the biological system, e.g., 
a synthetic genetic oscillator [H, |?], a coupled-oscillator 
system that represents cardiovascular phenomena [8] and 
cellular differentiation j9[. 

In the context of AD and OD the work by Koseska 
et al. pa deserves a special mention, because it shows 
that AD and OD differ in manifestation and origin. Most 
significantly, it reports a Turing type bifurcation that 
marks a direct continuous transition from AD to OD. 
The observations of Ref. [10] have later been verified 
in different systems and coupling schemes and a recent 
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burst of publications explore many aspects of the AD, 
OD and AD-OD transitions (see for example Refs. 0, 
Refs. 12], and Ref. 0). But, a continuous endeavor 
to find various aspects of AD and OD in networks (e.g., 
chimera death |14[), new systems (e.g., ecological system 
0), and new coupling schemes (e.g., amplitude depen¬ 
dent coupling 0 ) indicates the necessity and urgency 
of further research in this field. 

In the previous studies on oscillation suppression, all 
the variables of a coupled oscillators under study show 
either OD or AD, separately, for a certain parametric 
condition. Under no condition it happens that different 
variables of a same system show different oscillation sup¬ 
pression states 0. This is owing to the fact that either 
the inverse Hopf bifurcation (that leads to AD) or a sym¬ 
metry breaking bifurcation (that leads to OD) occur to 
all the variables of a system at a time. 

At this point we ask the following open question: Is 
it possible that different variables of a system (of coupled 
oscillators) show different types of oscillation suppres¬ 
sion states under an identical parametric condition ? In 
this paper we indeed identify this type of mixed mode 
oscillation suppression states in the paradigmatic Lorenz 
system and Lorenz-like system under several generic cou¬ 
pling configurations, which were earlier studied in the 
context of AD and OD. 

In the present paper we report two mixed mode os¬ 
cillation suppression states : (i) The mixed mode death 
(MMD) state, where, under a similar parametric condi¬ 
tion, a set of variables of a coupled oscillator system show 
OD, whereas the rest of the variables of the same cou¬ 
pled system show the AD state. In the chaotic Lorenz 
system under two different coupling schemes, namely the 
direct-indirect 00 and mean-field coupling 1, 0, 
we show that, while two of the variables (say, x and y) 
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show the OD state, the rest (i.e., the z variable) shows 
an AD state. Thus, unlike the AD or OD state, this 
MMD state is a variable selective mixture of stable IHSS 
(i.e., OD) and stable HSS (i.e, AD), (ii) In the second 
mixed death state, a nontrivial bistable AD (NAD) state 
(occurs in x and y variables) and a monostable AD state 
(occurs in the z variable) appear simultaneously-we refer 
this state as the MNAD state (i.e., mixture of the NAD 
and AD state). 

We identify that the symmetry of the system plays an 
important role behind the birth of the MMD and MNAD 
states: While the reflection symmetry breaking is respon¬ 
sible for the MNAD state, the breaking of both reflec¬ 
tion and translational symmetries are responsible for the 
MMD state. Through a rigorous bifurcation analysis, we 
establish the occurrence of the MMD and MNAD states 
and several transition routes associated with them. To 
establish the generality of our results we verify all the 
results in a Lorenz-like system, namely the chaotic Chen 
system [20]. To the best of our knowledge, existence of the 
MMD and MNAD states and the corresponding transi¬ 
tions associated with them have not been observed earlier 
for any other system or coupling configuration. Finally, 
we support our results through an electronic experiment 
that provides the first experimental evidence of the MMD 
and MNAD states. 

The rest of the paper is organized in the following 
manner: The next section considers the identical Lorenz 
systems under two different coupling schemes. Rigor¬ 
ous bifurcation analysis establish the occurrence of MMD 
and MNAD, and related transitions. Section inn reports 
the occurrence of mixed death states in Chen system, 
which is a Lorenz-like system. Experimental observa¬ 
tion of MMD and MNAD is reported in Sec. IIV1 Fi¬ 
nally, Sec. E| concludes the outcome and importance of 
the whole study. 

II. COUPLED LORENZ SYSTEMS 

A. Identical Lorenz systems interacting through 
direct-indirect coupling 

At first, we consider two identical chaotic Lorenz oscil¬ 
lators [2l[ interacting directly through diffusive coupling 
and indirectly through a common environment s, which is 
modeled as a damped dynamical system [13|, Qj| S, [231 . 
The mathematical model of the coupled system is given 
by 


Xl ,2 

= v(yi,2 - £1,2) + d(x 2,1 - x 1>2 ) + es, 

(la) 

m,2 

= (r- 21,2)2:1,2 - 2/1,2, 

(lb) 

ii ,2 

= £1,22/1,2 - bzi, 2> 

(lc) 

s 

_ ^ e(x 1 +x 2 ) 

(Id) 

: d is 

the diffusive coupling strength, and e 

is the en- 


vironmental coupling strength that controls the mutual 


interaction between the systems and environment, k is 
the damping factor of the environment (fc > 0) 1181. This 
coupling scheme was proposed by Resmi et al. [l8| as a 
general scheme to introduce AD in any coupled oscilla¬ 
tors. Although Ref. [l8| studied the response of a Lorenz 
system under this coupling scheme, but it could not iden¬ 
tify the OD state; only AD was reported there. Later the 
present authors [13| reported that this coupling scheme 
can also gives rise to OD and AD-OD transitions in non¬ 
linear oscillators. Thus, it is a generic coupling scheme to 
introduce AD and OD. Here, in the following, we inves¬ 
tigate the effect of this coupling scheme in inducing the 
mixed mode oscillation suppression states in the Lorenz 
system. 

Equation m has a trivial homogeneous steady state 
(HSS), which is the origin (0, 0, 0, 0, 0, 0, 0), and addition¬ 
ally two more coupling dependent nontrivial fixed points 
given by: 

(i ) Fmmss = (a: t ,2/ t ,-x\-y 1 ' ,^ t ,s t ), where = 

I / {(r— 1 )<T — 2 rfjfe + _ (<j+2rf)x t + _ (<7+2d)a: t2 f _ q 

Fmmss gives the mixed mode steady states (MMSS) as 
for these steady state we have inhomogeneity in x and y 
variables (i.e. ±x i = ^X 2 ,±yi = T 2 / 2 ) and homogeneity 
in the z variable (i.e., z\ = Z 2 ) • The stabilization of 
Tmmss results in the mixed mode death (MMD) state, 
because here OD occurs in x and y variables and AD 
occurs in the z variable. Note that Fmmss depends only 
upon d and independent of e and k. 

(i i) Fnhss = ( x*,y*,z*,x*,y*,z*,s *), where x* = 

1 / (fecr(r —1) —e 2 }b * _ (. ka+e 2 * )x * * _ (. ka+e 2 )x * 2 

Y /ccr+e 2 ’ U kcr 5 kerb ’ 

s* = Fnhss represents nontrivial homogeneous 

steady states (NHSS), stabilization of which gives rise to a 
novel nontrivial amplitude death (NAD) state (observed 
in x and y variables with ±x\ = ±£ 2 , ± 2/1 = ± 2 / 2 ), which 
is a nonzero bistable state and a monostable AD state 
(observed in the z variable with z\ = Z 2 )~we refer this 
state as mixed NAD and AD state, i.e., the MNAD state. 
It can be seen that, J~nhss depends upon e and k , and 
independent of d. 

As we notice here the symmetry of the system plays an 
important role behind the birth of Tmmss and J~nhss , 
and thus the emergence of the MMD and MNAD states. 
The uncoupled Lorenz oscillators, denoted by, say, C\ 
and £2 are invariant under the reflection about the z- 
axis and Eqs. m show the presence of a translational 
symmetry between the two Lorenz oscillators C\ and C 2 
under this coupling scheme. In the MNAD state, the re¬ 
lation between the oscillator variables clearly shows that 
the translational symmetry between C\ and C 2 is pre¬ 
served [Vi = V 2 {V = x, y ), z\ — Z 2 ] , however, the 
reflection symmetry about the z- axis is collapsed, as now 
Up 2 and —Hi ,2 are two different states. On the other 
hand, in the MMD state, the relation between the vari¬ 
ables of C\ and C 2 shows that V 2 = —Vi, which indicates 
the destruction of translational symmetry, and also, since 
now Vq 2 and — Vi ,2 represent two different states, we can 
say that the reflection symmetry about the 2 -axis is also 
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FIG. 1. (Color online) Bifurcation diagram with variation 
of d (using MATCONT) of two identical chaotic Lorenz sys¬ 
tems under direct-indirect coupling (e = 20, i.e., e > cpB 2 )• 
Grey (red) lines: Stable fixed points, Black lines: Unstable 
fixed points; open circle (blue): Unstable limit cycle, x and y 
variables show OD but the z variable shows AD (AD2). The 
MMD state (0D+AD2) (shown in shade) is created from ADI 
through a supercritical pitchfork bifurcation at PB1 with de¬ 
creasing d and destroyed at HBS3 (subcritical Hopf bifurca¬ 
tion point). Other parameters are k = 1, r = 28, a = 10, 

bUl 

broken. Thus, in the MMD state both the translational 
symmetry and the reflection symmetry are broken. 

Next, we theoretically analyze the stability of the sys¬ 
tem in order to explore the bifurcation scenarios. The 
characteristic equation of the system at the trivial HSS, 
(0,0,0,0,0,0,0), is given by, 

(6 + A) 2 F 2 (A)F 3 (A)=0, (2) 

where F 2 (A) = (A 2 +P^A+P^ 0 ), F|,(A) = (A 3 +P T2 A 2 + 
Pti\ + Pto) with P' TX = (l + cr + 2d), P^ 0 = ^(l — r) + 2d, 
Pt 2 = (1 + <j + fc), Pti = cr( 1 — r) + k( 1 + o) + e 2 , and 
Pro = e 2 + Ka(l — r). Since Eq.(|2]) is a seventh-order 
polynomial, it is difficult to derive the exact analytical 
expressions of all the eigenvalues. But we can predict 
the stability scenario of the trivial HSS from the proper¬ 
ties of the coefficients of the characteristic equation itself 
[24j . We find that all the eigenvalues at the trivial HSS 
have negative real part and thus give rise to ADI, when 
k > — (1 + cr), d > an d e > y/ka(r — 1). With 

decreasing strength of the coupling parameters (d and 
e) we observe different dynamical behaviors in the cou¬ 
pled system. A close inspection of the nontrivial fixed 
points reveals that Pmmss and Pnhss appear through 
a pitchfork bifurcation at dpBi and cpb 2 , respectively, 
where 


_ (r - 1)<7 
dpB 1 ^ ’ 

(3) 

eps2 = \Jko(r - 1 ). 

(4) 
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FIG. 2. (Color online) Bifurcation diagram of two identical 
chaotic Lorenz systems under direct-indirect coupling with 
variation of e. x and y show the NAD state and z shows an 
AD state (AD3). The MNAD state (shown in shade) emerges 
from ADI at PB2 (supercritical pitchfork bifurcation point) 
with decreasing e (d = 150, i.e., d > dps 1 ) and destroyed at 
HBS4 (subcritical Hopf bifurcation point). Other parameters 
are same as Fig. [I] 


To keep the uncoupled Lorenz systems in the chaotic re¬ 
gion, we set the system parameters at r = 28, b = | and 
cr = 10 throughout this paper. 

Figure [1] shows the one-dimensional bifurcation dia¬ 
gram with variable d using the MATCONT package [25| 
with an exemplary value k m 1 and e = 20 (note that 
e > €pb 2 = 16.432). It can be clearly seen that x\^ 
and 2/1 s 2 variables show the OD state [Figs. [1] (a) and[l] 
(b)], whereas, the zi } 2 variables show an AD state (AD2) 
[Fig. DJc)]; both these OD and AD states emerge from 
the amplitude death state (ADI) through a supercritical 
pitchfork bifurcation at dpsi = 135 [this value exactly 
matches with Eq. ©]• Thus, for e = 20, a decreasing 
d gives rise to an MMD state (OD+AD2) from an AD 
state (ADI). The MMD state becomes unstable through 
a subcritical Hopf bifurcation (HBS3). 

Next, we fix the value of d at d = 150 (d > dpBi), 
and vary e. Figure [2] shows that, with decreasing e, the 
bistable NAD state occurs in xi : 2 and 2 / 1,2 [Figs. El(a) and 
[2][b)] and an AD state (AD3) occurs in [Fig. EJc)] at 
tPB 2 = 16.432 [this value exactly matches with Eq.(j4])]; 
both these states emerge from the ADI state. Thus, a 
decreasing e (and a proper value of d) gives rise to an 
MNAD state (NAD+AD3) from the ADI state. The 
MNAD state is destroyed through a subcritical Hopf- 
bifurcation (HBS4). 

To identify several other bifurcation curves that mark 
the distinct regions of occurrence of oscillation suppres¬ 
sion states and their coexistence we consider the char¬ 
acteristic equation corresponding to the nontrivial fixed 
points (x l , y\ z l , Jx l , Jy l , z 1 , s l ) , where J = ± 1 and i = 
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FIG. 3. (Color online) Between HBS4 and PBS2 (a) OD and 
NAD and (b) AD2 and AD3 coexist, (a) Inset shows the 
zoomed in view of the NAD state, (c) NAD and OD and (d) 
AD3 and AD2 coexist between PBS2 and HBS1. HBS and 
PBS represent the subcritical Hopf and subcritical pitchfork 
bifurcation points, respectively. Other parameters are same 
as in Fig. [T] 


* or f, which is given by 

Fn(A)Fn(A) = 0, (5) 

where Fg,(A) = (A 3 +i$ ? A 2 +i$p+i$ 0 ), andF*(A) = 
(A 4 + P^ 3 A 3 + P^ 2 A 2 + P^iA + P 1 N q ) with P^ 2 = & -\~2d-\- 
b + 1, Pjyi — b + (b + 1)(<j + 2d) + — r), 0 = 

(cr + 2d)(b + x z2 ) + a{x l y l + b(z l -r)}, P^ 3 = fc + 6 + l + cr, 
P l N2 = fc(6 + 1 + cr) + cr{b + 1 - r + £ z ) + b + a;* 2 + e 2 , 
P^i = fc{cr(6 + 1 — r + z l ) + b + x t2 } + a {6(1 — r + 
2 d) + x l [x l + y 1 )} + e 2 (6 + 1) and P^ 0 = e 2 (6 + ad 2 ) + 
ka{b( 1 — r + 2 d) + ad (ad + ?/*)}. From Eq. (j5j) with J = 1 
and i = * (i.e., for the fixed point Pnhss ) we find that 
d appears only in the term F^(A), i.e., it controls only 
three eigenvalues. Similarly from Eq. © with J = — 1 
and i = f (i.e., for the fixed point Pm mss) we conclude 
that behavior of the four eigenvalues are controlled by e 
and k as they appear only in the term *n(A)- 

For d < dpBii k = 1 and €pbs 2 > ^hbsai the MNAD 
state appears through a pitchfork bifurcation at PBS2. 
The analytical expression is obtained by putting Pj^ Q = 0 
0 and is given by 

dFBS2 = -2(6 +*> 2 )-■ <6) 

Another Hopf bifurcation point HBS2 is observed for 
lower values of d. To derive the locus we set (P^ 1 P^ 2 — 
Cvo) =00 and get 


dHBS2 = 


-P* + a/p* 2 - 4A*C* 


(7) 



FIG. 4. (Color online) Two-parameter bifurcation diagram 
(using XPPAUT) in e — d space. Other parameters: k = 1, 
r = 28, cr = 10 and b = |. 


where, A* — 4(6T 1), B * — 2(6Tl)(2crT 6 T 1) T 2<t(z* — 
r), C* = <j 2 (b + 1 - r + z*) + a{(6 + l) 2 - r + - 

x*y*} + {b+x* 2 )(b + l ). When d < dp b, depending upon 
the choice of k and d values, the MNAD state appears 
either through PBS2 or HBS2. Figures [3fa, b) show that 
for d = 60, i.e., d < dpsi and k = 1, the ADI state 
disappears; coexistence of NAD (AD3) and OD (AD2) 
is observed between PBS2 and HBS4. For e < ep# 2 , 
the ADI state vanishes and we get two possible routes 
to the MMD state: (i) Pitchfork bifurcation (PBS1) and 
(ii) Hopf bifurcation (HBS1). With the proper choice of 
k and e we can select one of these routes to the MMD 
state. Figures [3{c, d) show the Hopf bifurcation route 
to MMD for e = 15, i.e., e < cpb 2 and k = 1. To get 
the exact locus of PBS1 and HBS1 we set P^ 0 = 0 and 

( P N 1 P N 2 P N 3 - pi Ni - P N0 P N3 2 ) = °> respectively, and 
get 


tPBS 1 = 


£HBS 1 = 


I ka{b(r — 1 — A) — ad (#t + yt)} 


b + ad 


'-Pt + \/P t2 - 4AtCt 
2At : 


(8a) 

(8b) 


where, A* = (b + 1 )(k + cr), P' = (b • l)[ptjv 3 (ptjv 2 - 
e 2 ) - 2{Pt jvl - e 2 (b + 1)}] + P* N3 [{Pl N i - e 2 (6 + 1)} - 
PW& + ^ 2 )], Ct = Pt JV3 [(P t iV2 - e 2 ){-Ptjvi - e 2 (& + 
1)} - PtjvajFtivo - e 2 {b + A)}] - {P^ N1 -e 2 (b + 1)} 2 . 

The loci of HBS4 and HBS3 could not be found in 
the closed form; thus, to present a complete bifurcation 
scenario we compute the two parameter bifurcation dia¬ 
gram (Fig. [4]) in the e — d space with k = 1 [i.e., with 
k > — (1 + cr)] using the XPPAUT package (26|, which 
exactly agree with our theoretically obtained bifurcation 


2A : 


curves. 
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B. Identical Lorenz systems interacting through 
mean-field diffusive coupling 

To verify that the mixed mode oscillation suppression 
states are not limited to the direct-indirect coupling only, 
we consider another generic coupling scheme in the con¬ 
text of AD and OD, namely the mean-field diffusive cou¬ 
pling El, and investigate the occurrence of MMD and 
MNAD states. We consider the following two identical 
mean-field coupled Lorenz systems: 


xi ,2 = cr(y i j2 - £ 1 , 2 ) A £ 


Q(x i + x 2 ) 


~ £l,2 


2/1,2 = (r- Zi t2 )xi ,2 ~ 2/1,2, 
Zl ,2 = £ 1 , 22 / 1,2 - bz 1 , 2 . 


(9a) 

(9b) 

(9c) 


Here e is the coupling strength and the control pa¬ 
rameter Q determines the density of mean-field (0 < 
Q < 1) [27[. From Eqs. (j9j) we can see that the origin 
(0, 0, 0, 0, 0, 0) is the homogeneous steady state (HSS). 
Also, we have two more coupling-dependent nontriv¬ 


ial fixed points: (i) Fmmss = (x\y^z\—x^ 
where = =b 

x^ 2 (a +e) 


-y' 


,g, 


b(ar-a-e) 
cr+e 5 


2/t = zt = 


bcr 

X* = =b 


(ii) Fnhss = (x*. 

y * 


X 


b{a(r-l)-e(l-Q)} 
cr +e(l —Q) 

{cr+e(l — Q)} 


y*, z*,x*z*), where 

_ x*{a+e(l — Q)} 


Z~ = 


bcr 


The nontrivial fixed points Fmmss and Fnhss 
emerge due to the symmetry breaking pitchfork bifur¬ 
cations at 6 pb 2 and ep#i, respectively, 


£pb l = 


j(r — 1) 

1 -Q ’ 
£pb2 = &(r — 1 ). 


(10a) 

(10b) 


To explore the complete bifurcation scenario, we write 
the characteristic equation of the system at the nontrivial 
fixed points (x % , y l , z l , Jx l , Jy l , z l ) , where J = ±1 and 
i = * or f, as 


F£(A)F" 3 (A) = 0 


( 11 ) 


where, F$(A) = (A 3 + U% 2 \ 2 A U^X + U% 0 ) with 
l = 'or ", Uff 2 = b A 1 + cr + e(l — Q), U ^ = 
b A (6 + l)(cr + e - eQ) + x l 2 - <r(r - 2 *), {/$- 0 = 
cra:* 2 /*-&cr(r- 2 :*) + (cr+e-e(5)(a;* 2 +6), U'^ 2 = a+e+b+l, 
U'ni = P + e)(b + 1 ) + b + x i2 + a(z i - r), £/$ 0 = 
(cr + e)(6 + a:* 2 ) + ax l y l A ba(z l — r). 

Using a similar approach adopted in the Sec. ESI we 
derive the locus of the bifurcation curves and they are 
given by 


£pbs = 


Qhbs i = 


y[rQ' — 4 + \/r 2 Q ' 2 + 8Qr] 


-j?|- y/B?-AA{c\ 

2 A\ 


(12a) 

(12b) 



FIG. 5. (Color online) Bifurcation diagram of two identical 
chaotic Lorenz attractors coupled through mean-field diffu¬ 
sion at Q — 0.5. Here MNAD (NAD+AD3) is created at 
PB1; The MMD state (OD+AD2) created at PBS is always 
accompanied by the MNAD state; PBS: Subcritical pitchfork 
bifurcation, (a) In x variable: ADI to NAD transition (at 
PB1) and co-existence of OD and NAD (between PBS and 
HBS1). Same is for the y variable (not shown here), (b) In z 
variable: ADI to AD3 transition (at PB1) and co-existence of 
AD2 and AD3 (between PBS and HBS1) is observed. (r = 28, 
<j — 10, b = |). 


Here Q' = (2 -Q), A\=e 2 {b + 1), B] = -(<x + e)[^I + 
2 c(b “h 1)] A ore — {b A l)^e, = [6~hl — cr — c A 

( 1+<7+ ^( ff+6 ) ]£ t2 -crr(l+o’+e)+6(6+l) + (6+l) 2 (c r +e) + 
(b + l)(cr + e) 2 . Here cpbs gives the coupling strength at 
which the MMD state emerges (due to the stabilization 
of Fmmss) through a subcritical pitchfork bifurcation. 
From Eq. (jlQb) ) and Eq. (j!2a|) it is clear that although 
the emerging point of Fmmss is independent of Q, but 
its stabilization, i.e., the creation of the MMD state, is 
controlled by Q. 

Figures [5fa) and[5^b) show the bifurcation diagram of 
xi ? 2 and 2 i } 2 , respectively for Q = 0.5 [ 2 / 1,2 behaves in 
a similar way as ^ 1,2 and thus not shown in the figure]. 
With this coupling scheme we obtain all the oscillation 
quenching states, namely the MMD (i.e., OD + AD2) 
and the MNAD (i.e., NAD + AD3) state. It is notewor¬ 
thy that the MMD state is always accompanied by the 
MNAD state. Also, here the ADI to MMD transition 
is not possible, as for any Q > 0 one has cpbi > £pb 2 - 
However, the direct transition from ADI to MNAD takes 
place at cpbi- 

The complete bifurcation scenario is shown in Fig.[6]in 
the e — Q space. Here the intersection of HBS1 curve with 
PBS organizes the coexisting MMD and MNAD state, 
and the region bounded by HBS2 and PB1 curve orga¬ 
nizes the occurrence of the MNAD state. The horizontal 
dotted line indicates the density of the mean field for 
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FIG. 6. (Color online) Two-parameter bifurcation diagram of 
two identical chaotic Lorenz systems coupled through mean- 
field diffusion. The horizontal dotted line indicates the Q 
value for which Fig. [5] is drawn (r = 28, a = 10, b = |). 


which Fig. [5] is drawn. 


III. CHEN SYSTEM 


A. Identical Chen systems interacting through 
direct-indirect coupling 


Next, we verify the generality of the occurrence of the 
mixed mode oscillation suppression states in a chaotic 
Lorenz-like system, namely the Chen system 120]. Math¬ 
ematical model of two identical Chen systems under the 
direct-indirect coupling scheme is given by 


*1,2 = a(yi,2 - an,2) + d(x 2,t - an,2) + es, 
2/1,2 = (c - a)an, 2 - an, 2 zi,2 + ct/i,2, 

*1,2 = 2:1,22/1,2 - &Zl,2, 
s=-ks- e(x '+ X2) . 


(13a) 

(13b) 

(13c) 

(13d) 


Here a > 0, c ( 2 c > a) and b > 0 are the system pa¬ 
rameters. In addition to the trivial homogeneous steady 
state (HSS), i.e., the origin ( 0 , 0 , 0 , 0 , 0 , 0 , 0 ), the system 
has two more coupling dependent nontrivial fixed points 
(i) Fmmss = {x \y\z\-x\-y\z\s^), where x f = 

1 / b{2c(q+d)-a I 2 } + _ a; t (q+2d) + _ x^(a+2d) f _ q 

and (ii) Fnhss = (,x*,y*,z*,x*,y*,z*,s*), where x* = 

I j b{ka(c—a)+c(ka+e 2 )} * x* (ka-\-e 2 ) * x* 2 (fca+e 2 ) 

Y ka-\-e 2 7 y ka 7 kab 7 

s* = A close inspection of the nontrivial fixed 

points reveals that 


dps 


a 2 — 2ac 


2c 


(14) 


where dpB gives the coupling strength at which J~mmss 
emerges. 

From the detailed analysis of the characteristic equa¬ 
tion we get the stability condition of the trivial HSS as 
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FIG. 7. (Color online) Bifurcation diagram of two identical 
chaotic Chen systems coupled through direct-indirect cou¬ 
pling. (a,b) With e = 5, the MMD state (OD+AD2) born 
through subcritical Hopf bifurcation (HBS1). NAD (AD3) 
state exists between HB5 and HB4. (a) Inset shows the 
zoomed in view of the OD state. (c,d) With d = 4 MNAD 
is noticed, (c) NAD and (d) AD3 state appear at HB4. PB, 
PBS, HB and HBS denote the pitchfork, subcritical pitchfork, 
Hopf, subcritical Hopf bifurcation points, respectively. Other 
parameters are k = 1, a = 40, c = 28, b = 3. 


k > (c — a), e 2 < fe ( Q ~ 2ac ) an q ^ < q2 ~^ qc . Consider¬ 
ing the conditions for the system parameters (i.e. a > 0, 
2 c > a and b > 0) of the chaotic Chen attractor we get 
(a 2 — 2 ac) < 0 for all possible set of parameter values. So 
for any positive &, the stability condition for the trivial 
HSS is satisfied when e 2 < 0, i.e., for imaginary e val¬ 
ues. These conditions clearly show that for any positive 
fc, the trivial HSS remains always unstable and the AD 
state, that arises due to stabilization of the trivial HSS, 
never appears. The detailed bifurcation scenario of the 
system is shown in Fig. [7] with an exemplary value k = 1. 
Figures [3(a) and[7fb) show the bifurcation structure for 
e = 5; here we can see the presence of MMD (OD+AD2) 
and MNAD (NAD+AD3). To show the complete bifur¬ 
cation structure we also consider the negative d values 
(and later, also negative e values). Figures [7] (c) and [7| 
(d) show the bifurcation for varying e and fixed d = 4; it 
shows the presence of MNAD (NAD+AD3) state in the 
coupled system. 


B. Identical Chen systems interacting through 
mean-field diffusive coupling 


Next, we verify the occurrence of MMD and MNAD 
states in two identical Chen systems under the mean- 
field diffusive coupling scheme. The mathematical model 
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FIG. 8. (Color online) Bifurcation diagram of two identical 
chaotic Chen systems coupled through mean-field diffusion. 
The stabilization of Rmmss (Fnhss) gives rise to (a) OD 
(NAD) and (b) AD2 (AD3) state. Other parameters are Q = 
0.05, a = 40, c = 28, b — 3. 


of the coupled system is given by 


, , . (Q(x 1 +X 2 ) \ , . 

X\,2 = a{yi,2 - a:i,2) + e| --- x\ J, ( 15 a) 

yi,2 = (c - a)x i >2 - xi,2Zi,2 + cy i j2 , ( 15 b) 

Zl,2 = Xl,2Vl,2 - bzi t 2- ( 15 c) 

The Eq. (ITK 1 ) has the trivial fixed point ( 0 , 0 , 0 , 0 , 0 , 0 ) 
and two more coupling dependent nontrivial fixed 
points (i) Tmmss = (x f , z f , -x f , -y\z*) where 

f _ j_ / ab(c—a)+cb(a+e)~ + _ xUji+e) f _ 

Y a+e i y a 5 

- lfc +£) a nd (ii) Tnhss_ = (x*,y*,z*,x*,y*,z*), where 

* _ i / b{2ac-a 2 +ce{l-Q)} * _ (a+e(l-Q)} * _ 

x ~ a+e(l-Q) ’ y — a ’ Z “ 

{a+ aft (1 ~ Q)} - Fmmss and Fnhss born through the 
pitchfork bifurcation at 6pb2 and ep#i, respectively, 
where 


a 2 — 2ac 
c(l - Q)’ 


(16) 


a 2 — 2 ac 

tPB2 = - 

C 


(17) 


The stabilization of Tmmss and Tnhss gives rise to 
the MMD and MNAD states, respectively. Figure [Sfa) 
and[8fb) show the bifurcation diagram of and zi^, 
respectively for Q = 0.05. Both the MMD (OD+AD2) 
and MNAD (NAD+AD3) states appear through subcrit- 
ical hopf bifurcation at HBS2 and HBS4, respectively. 


IV. EXPERIMENT 

We experimentally verify the occurrence of MMD and 
MNAD states in the identical Lorenz attractor inter¬ 



FIG. 9. (Color online) Experimental circuit diagram of two 
coupled Lorenz systems under direct-indirect coupling. Al- 
A8, C1-C6 are realized with TL082 op-amps. All the unla¬ 
beled resistors have value Ri — 100 kH, C=10 nF, V a — 3 
volt, R — 10 kQ, Rb — 39 kQ, R r — 1.07 kfl, R a — 10 kfl, 
Rk = 100 kQ. ±15 volt power supplies are used; resistors (ca¬ 
pacitors) have ±5% (±1%) tolerance. Box denoted by “B” are 
op-amp based buffers; inverters are realized with the unity- 
gain non-inverting op-amps. Insets (in the middle part) shows 
the experimental attractors from the uncoupled Lorenz oscil¬ 
lators Ci (left) and £2 (right) (y- axis: 1 v/div, x-axis: 1 
v/div). 


acting through direct-indirect coupling. To implement 
the practical electronic circuit we have rescaled [28| the 

variables of Eq. JD using x* = (y7 ) x i » Vi = 

Zi = (^)zi, Si = (y7) s * where i = 1,2. Then the modi¬ 
fied equations become 


xi ,2 = <?(yi,2 ~ xi,2) + d{x 2,1 - xi >2 ) + es, 
r 

y 1,2 = g(3 — Zi t2 )xi t 2 - y 1,2) 
z 1,2 = Xl, 2 yi,2 - bz 1 , 2 , 

£^+£5) 

2 


(18a) 

(18b) 

(18c) 

(18d) 


Figure. [9] represents the electronic circuit of the cou¬ 
pled Lorenz systems [28] with the direct-indirect coupling 
given by Eqs. (fl8j) . A1-A4 and A5-A8 are used to realize 
the individual Lorenz oscillators C\ and £ 2 , respectively. 
The subunit realized with op-amps C5-C6 produces the 
diffusive coupling part, and the subunit consists of op- 
amps C1-C4 mimics the environmental coupling part. 
The voltage equation of the circuit of Fig. [9] can be writ¬ 
ten as follows: 
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FIG. 10. (Color online) Transition from limit cycle (LC) 
to MMD: (a, c, e) Experimental real time traces of VA, V y i 
and Vzi along with the (b, d, f) numerical time series plots 
of Xi, y% and z% where i — 1,2. R e — 1.127 kQ (e = 11.27): 
Left panels show LC for Rd = 4.55 kQ (d — 45.5). Right 
panels show MMD for Rd — 145.6 Q (d = 1.456). Both 
the experimental and numerical results clearly show that the 
MMD state is a mixed state of OD (observed in x, y and V x , 
V y ) and AD (AD2) state (observed in z and 14). For other 
parameters see text. Left panel, y axis: (a) 100 mv/div (c) 
340 mv/div (e) 128 mv/div; x axis: 800 /xs/div. Right panel, 
y axis: (a) 100 mv/div (c) 300 mv/div (e) 128 mv/div; x axis: 
10 ms/div. 


10 5 

10 5 CV xi = —(V yi -V xi ) + Xc, 

lX(j 

(19a) 

10 4 

10 5 CV yi = — (V a - V zi )V xi - Vyi, 

ilf 

(19b) 

10 5 

10 5 CV zi = v xi v yi - —v zi , 

Kb 

(19c) 

10 5 cv s = 105 Vs Re ( Vxl + Vx2 ). 

R k 100 V 2 ) 

(19d) 

Where I c = $ fr(V xj - V xi ) + j^V s with i, 

3 = 1, 2 


and i 4 j. Equation m is normalized with respect to 
10 5 C and thus now becomes equivalent to Eq. (fl8l) for the 
following normalized parameters: u = r = 1Q | C , e = 


r l = RsL k= 10- fT= 1 = 10- L 

100 ’ a 100 ’ ^ Rk’ Ra ’ 3 Rr ’ 

ry» . jTc i n i . ^7 ]_ ry . Z i Q T~l cl C Vs 

dj i v + ’ y* 1 v /C ' % v # 5 * v + 

v sat -r* sat i t-» v sat rv sat 

resistances Rd, R e and Rk control the diffusive coupling 


iQl v ~ 3 
-. Thus, the 
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FIG. 11. (Color online) Transition from limit cycle (LC) 
to MNAD: (a, c, e) Experimental real time traces of VA, Vyi 
and Vzi along with the (b, d, f) numerical time series plots 
of Xi, yi and zi where i — 1,2. Rd — 8.5 kQ (d m 85): Left 
panels show LC for R e — 0.84 kQ (e = 8.4). Right panels 
show MNAD for R e = 1.4 kH (e — 14). For other parameters 
see text. Left panel, y axis: (a) 156 mv/div (c) 300 mv/div 
(e) 128 mv/div; x axis: 780 /xs/div. Right panel, y axis: (a) 
50 mv/div (c) 350 mv/div (e) 100 mv/div; x axis: 10 ms/div. 


strength (d), environmental coupling strength (e) and the 
damping factor of the environment (k), respectively. V sa t 
is the op-amp saturation voltage. We choose C=10 nF, 
R r = 1.07 ktt (r = 28), R a = 10 kQ (a = 10), R b = 
39 kVl (b = 2.6), and V a = 3 volt. These particular 
choice of parameter values make the system represented 
by Eq. (fl9| ) equivalent to that given by Eq. (fl8l) and keep 
the uncoupled Lorenz systems in the chaotic region. 

At first, we take R e = 1.127 kQ (e = 11.27), Rk = 
100 k Vl (fc = 1) and observe a continuous transition 
from the limit cycle to MMD for decreasing Rd (d). In 
Figs. flOl a). (c), (e) using the experimental snapshots 
of the wave forms [taken using a digital storage oscillo¬ 
scope (Agilent, DSO-X 2024A, 200MHz, 2 GSa/s)], we 
experimentally observed two distinct dynamical regions 
for two different Rd (d) values. (1) Limit cycle (LC): 
For Rd = 4.55 kfl i.e. d = 45.5 [Figs. HQ](a), (c), (e) 
left panel]. (2) MMD: For Rd = 145.6 Vt i.e., d = 1.456 
[Figs. flOl a). (c), (e) right panel]. We define this MMD 
state (see Sec. EU as the mixture of OD and AD state. 
The experimental real time traces clearly show this si¬ 
multaneous occurrence of OD in V x [Fig. [TO] (a) right 
panel] and V y [Fig. [TOlT c) right panel] and AD (AD2) 
state in V z [Fig. [TOf e) right panel]. The numerical time 
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series plots (using the fourth-order Runge-Kutta method 
with 0.01 step size) for the equivalent parameter values 
are shown in Figs. flOfb).(d) and (f), which clearly shows 
the qualitative agreement between the experimental and 
numerical results. However, the slight mismatch between 
the experimental and numerical results may be due to 
the possible parameter mismatch and fluctuations that 
are inevitable in a real circuit. 

Next, we set Rd = 8.5 kD (d = 85), Rk = 100 kD 
(k = 1) and observe a continuous transition from the 
limit cycle to MNAD for increasing R e . The experimental 
results are shown in Figs. fTTla). (c) and (e). Here the 
following observations are made (1) LC: For R e = 0.84 kD 
(e = 8.4) [Figs. DjJa), (c), ( e ) left panel]. (2) MNAD: 
For R e = 1.4 kfl (e = 14) [Fig. QTJa), (c), (e), right 
panel]. The numerical time series plots are demonstrated 
in Figs.[l]Jb),(d),(f) for the equivalent parameter values. 
Here also, Fig. [lT]clearly shows that the MNAD state, as 
defined in Sec. un is a mixed state of a bistable NAD and 
AD state. In the experimental study the bistability of 
the NAD state is found by a random parameter sweeping 
around R e = 1.4 kfl, and the same is verified in numerical 
simulations using the proper initial conditions (not shown 
in Fig. fill) . 


V. CONCLUSION 

We have reported a new cooperative dynamical state, 
namely the mixed mode oscillation suppression state , 
where different set of variables of a system of coupled 
oscillators show different types of oscillation suppression 
states under the same parametric condition. We identify 
two types of this state in coupled chaotic Lorenz oscilla¬ 
tors: One is called the mixed mode death (MMD) state, 


where OD and AD occurs simultaneously to different set 
of variables, and the other one (called the MNAD state) 
is the mixed variable selective state of nontrivial bistable 
AD and a monostable AD. To show the generality of the 
results we consider two generic coupling schemes, namely 
direct-indirect coupling and mean-field coupling, which 
were studied earlier in the context of AD and OD. Also, 
we verify the results in the coupled Chen system, which 
is a Lorenz-like system. Through rigorous bifurcation 
analyses we find all the transition routes to these mixed 
oscillation suppression states and map them in parame¬ 
ter space. We identify the underlying symmetry breaking 
that leads to the MMD and MNAD states. Finally, we 
report the first experimental observation of the MMD 
and MNAD state using an electronic circuit experiment. 

The present study may have applications in many real 
systems, such as laser (29[ and geomagnetic [30] systems, 
whose models mimic the Lorenz system (under some 
proper transformations). Take for example of the laser 
system modeled by Maxwell-Bloch equation [2j| where 
the variables related to the electric-field and polarization 
mimic x and y variables, respectively, whereas the vari¬ 
able related to population inversion mimics the z variable 
of the Lorenz system. Thus, we believe that the results of 
the present study can be extended to other “Lorenz-like” 
physical systems and may be useful in understanding of 
those systems. 
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